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PROBLEM 


Predict  the  mean  concentration  distribution  of  polymer  far  downstream  of  the 
injection  of  a  polymer  solution  into  a  two-dimensional,  constant-pressure  turbulent 
boundary  layer  on  a  flat  plate. 


RESULTS 

The  available  experimental  data  and  empirical  relations  for  the  mixing  of  a  pas¬ 
sive  scalar  quantity,  such  as  smoke  or  dye  after  injection  into  a  turbulent  boundary 
layer,  appear  to  be  consistent  with  the  outer-layer  similarity  relations  for  high 
Reynolds  numbers.  The  final-zone  similarity  laws  of  mixing-which  relate  mean  con¬ 
centrations,  free-stream  speed,  boundary-layer  thickness,  and  wall  shear  stress— are 
developed  within  the  framework  of  equilibrium  boundary-layer  flow,  assuming  the 
velocity-defect  law  of  Coles  and  using  some  experimental  results  of  Poreh  and 
Cermak,  Fiedler  and  Head,  and  Kibens.  The  analogy  of  polymeric  friction  reduction 
to  a  “negative  roughness”  implies  that  the  polymer  can  be  considered  to  be  a  passive 
scalar  outside  the  thin  inner  layer,  so  that  the  major  mixing  process  is  unchanged. 
The  analogy  is  expected  to  be  valid  in  some  practical  ranges  of  concentration  and 
wall  shear  stress  whose  limits  must  be  determined  experimentally. 

The  preliminary  results  of  mixing  experiments  in  an  open-channel  boundary 
layer,  for  injection  of  dyed  water  or  dyed  polymer  solution,  appear  to  be  consistent 
with  the  predictions  based  on  the  negative  roughness  analogy. 


RECOMMENDATIONS 

The  final-zone  similarity  relations  for  mixing  developed  here  should  be  tested 
by  further  experiments  in  nearly  constant-pressure  flows,  with  attention  to  truly 
two-dimensional  conditions.  The  goal  should  be  to  determine  the  ranges  of  polymer 
concentration,  wall  shear  stresses,  and  slot  distances  for  which  the  relations  are 
valid.  The  theoretical  predictions  of  mixing  and  friction  reduction  should  be  com¬ 
bined  to  predict  polymer  supply  rates  and  friction  reduction  for  injection  systems. 

In  comparing  such  predictions  with  experiments,  attention  must  be  given  to  the 
possible  importance  of  the  earlier  mixing  zones  and  the  concentration  fluctuations 
about  the  mean,  which  are  not  considered  in  this  work. 
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NOMENCLATURE 


B  Constant  in  logarithmic  law  of  the  wall:  B  -  B0  +  A B 
Bq  Smooth-wall  value  of  B  in  Newtonian  fluids 
c  Polymer  concentration  in  ppm  by  weight 
Cf  2rw/pU 2 

Cj  Value  of  c  in  injection-slot  fluid 

Cq  Volume-flux  average  value  of  c  in  the  boundary  layer:  QjCj/Q 

ct  Volume-flux  average  value  of  c  in  the  turbulent  fluid:  QjCf/Q, 
c  Time  average  value  of  c  at  a  point 
c 1  Turbulent  conditional  average  value  of  c  at  a  point 
cw  Wall  value  of  c 

G  Equilibrium  profile  parameter:  fQ  [(£/ -  u)/uT]2  drj/f0  ( U-u)/uTdi ? 
G  j  Clauser  wake  function 

H  'Profile  shape  factor:  //  ^  8*/ 6 

I  'Detector  function:  1  in  turbulent  flow;  0  in  irrotationa'l  flow 
k  Von  Karman  “universal”  constant;  k  «  0.41 
m  'Polymer  solution  mass  injection  rate  per  unit  span 
Q  Volume  flux  of  boundary-layer  fluid  per  unit  span 
Qj  Volume  flux  of  polymer  solution  per  unit  span 
Qt  Volume  flux  of  turbulent  fluid  per  unit  span 
t  Time 

u  x-component  of  the  velocity  vector 
uT  Local  value  of  the  friction  velocity:  y/rw/p 
u  Time  average  value  of  u  at  a  point 
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1 

If 

:  j 

id  Turbulent  conditional  average  value  of  u  at  a  point 
U  Free-stream  speed 
w(ri)  Coles’  wake  function 

x,y,  z  Rectangular  coordinates;  x  measured  along  surface  in  flow  direction ;y  nor¬ 
mal  distance  from  surface;  z  spanwise  position 

Y  Parameter  in  the  y  distribution 

7  Turbulence  internrittency  function:  / 

5  Boundary-layer  thickness  at  which  u  -  0.99 U 
5C1  Boundary-layer  thickness  used  by  Clauser 
8n  Boundary-layer  thickness  in  Coles’  law  of  the  wake 
5*  Boundary-layer  displacement  thickness:  (1 IU)  fQ  (U-u)dy 
A B  Change  in  B  with  wall  roughness  or  friction  reduction 

V  y/&„ 

0  Boundary-layer  momentum  thickness:  (\/(J)2  fQ  u(U-u)dy 
X  Value  of  y  at  which  c/cw  =  0.5 
v  Fluid  kinematic  viscosity 
n  Profile  parameter  in  Coles’  law  of  the  wake 
p  Water  (or  polymer  solution)  density,  a  constant 
o  Parameter  in  the  y  distribution 
tw  Local  wall  shear  stress 

Note:  In  the  cases  of  c,  cw ,  u,  uc ,  I ,  In,  and  Iiic,  the  overbar  indicates 
the  time-average  value  at  a  point  in  space  ( cr ,  Y ,  and  id  are  exceptions, 
which  are  explained  above). 
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INTRODUCTION 


The  hydrodynamic  frictional  drag  of  a  marine  vessel,  as  is  well  known,  can  be 
reduced  by  injecting  a  solution  of  any  of  various  long-chain  polymers  into  the  bound¬ 
ary  layer  on  the  vessel’s  hull  (Ref.  1  and  2).  The  polymer  solution  is  greatly  diluted 
by  turbulent  mixing  in  the  boundary  layer.  The  resultant  mean  concentrations  near 
the  hull  surface  along  most  of  the  vessel’s  length  should  be  close  to  certain  optimum 
values  if  the  maximum  friction  reduction  is  to  be  obtained  for  the  minimum  ex¬ 
penditure  of  polymer.  Thus,  the  prediction  of  boundary-layer  mixing  becomes  im¬ 
portant  to  the  evaluation  of  the  feasibility  of  drag-reduction  applications. 

Until  recently,  estimates  of  dilution  were  made  with  what  might  be  called  a 
perfect-mixing  model.  The  gradient  in  mean  concentration  across  the  boundary 
layer  was  neglected,  and  the  local  uniform  concentration  at  each  station  was  assigned 
according  to  the  boundary-layer  volume  flux  at  that  station  and  the  polymer  mass 
conservation  equation.  In  the  earliest  work,  that  volume  flux  was  estimated  simply 
by  integrating  the  mean  velocity  profile  from  the  wall  out  to  a  conventional  defini¬ 
tion  of  boundary-layer  thickness;  any  changes  in  velocity  profile  and  thickness  due 
to  the  polymer  solution  were  ignored.  Later,  the  method  of  estimation  was  improved 
by  considering  the  reduction  of  boundary-layer  thickness.  In  the  present  work,  major 
improvements  are  made  by  considering  the  intermittency  of  turbulence  and  polymer 
concentration,  the  profile  of  mean  concentration,  and  the  polymer  effects  on  the 
boundary-layer  mean  velocity  profile. 

Figure  1,  reproduced  from  Ref.  3,  illustrates  the  intermittency  of  a  turbulent 
boundary  layer  in  its  outer  part.  The  smoke-filled  boundary  layer  was  illuminated 
through  a  narrow  slit.  Fiedler  and  Head  (Ref.  3)  concluded  that  “smoke-filled  and 
smoke-free  regions  can  be  identified  with  regions  of  turbulent  and  non-turbulent 
flow.”  This  intermittency  must  be  considered  in  a  realistic  model  of  boundary-layer 
mixing. 

Real  understanding  of  the  mechanics  of  turbulent  boundary  layers  is  so  frag¬ 
mentary  that  idealized  models  of  behavior  must  be  relied  on.  For  example,  “eddy 
viscosity”  models  liken  the  turbulent  momentum  transfer  by  random  eddy  motion 
to  the  viscous  transfer  in  laminar  flow  by  molecular  agitation.  “Similarity”  models 
reduce  the  number  of  variables  that  have  to  be  considered.  Usually,  such  models 
are  accurate  only  within  the  range  of  the  experiments  from  which  the  similarity  laws 
are  deduced.  In  the  case  of  turbulent  boundary  layers  in  Newtonian  fluids,  a  well- 
developed  similarity  analysis  is  based  on  the  two-layer  concept  (Ref.  4).  For  large 
boundary-layer  thickness  Reynolds  numbers,  such  as  those  that  occur  on  marine 


FIG.  1.  Smoke  in  a  Turbulent  Boundary  Layer  (from  Ref.  3). 


vessels,  the  two-layer  model  is  particularly  useful  since  the  velocity-defect  similarity 
law  of  the  outer  layer  can  be  used  across  the  entire  boundary  layer  with  a  good  de¬ 
gree  of  accuracy. 

It  is  particularly  simple  to  incorporate  the  effects  of  friction  reduction  into 
boundary-layer  development  calculations  if  it  is  assumed  that  the  important  polymer 
effects  are  confined  to  the  inner  layer  (the  so-called  negative  roughness  analogy). 

This  assumption  should  be  valid  for  some  practical  range  of  operation,  and  several 
analyses  of  turbulent  boundary-layer  development  in  polymer  solutions  have  been 
based  on  it  (Ref.  2,  5  and  6).  In  the  present  work,  this  approach  is  extended  to  ob¬ 
tain  the  similarity  relation  for  the  two-dimensional,  far-downstream  mixing  of  any 
passive  scalar  quantity  injected  into  a  boundary  layer  with  friction  reduction.  Since 
that  relation  is  an  outer-layer  law,  it  is  reasonable  to  expect  that  it  also  applies  to  the 
far-downstream  mixing  of  injected  polymer  solution,  in  accordance  with  the  negative 
roughness  analogy. 

The  interest  in  predicting  the  boundary-layer  mixing  is  a  very  practical  one, 
since  the  effectiveness  of  a  polymer  solution  in  reducing  friction  depends  mainly 
upon  its  concentration  and  the  wall  shear  stress.  However,  experience  suggests  that 
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the  local  friction-reduction  effectiveness  of  a  boundary-layer  diluted  solution  of  a 
given  mean  concentration  may  be  much  less  than  that  of  a  homogeneous  solution  of 
the  same  concentration.  Such  behavior  may  be  expected  if  most  of  the  polymer  con¬ 
tributing  to  the  local  mean  concentration  is  in  the  form  of  blobs  of  small  volume- 
fraction  but  very  high  concentration. 

Thus,  the  limitations  of  the  present  results  must  be  borne  in  mind  when  they 
are  combined  with  boundary-layer  theory  to  predict  the  performance  of  injection 
systems.  In  particular,  the  present  work  predicts  only  the  far-downstream  time- 
average  concentrations  and  does  not  predict  the  intensity  of  the  concentration  fluc¬ 
tuations  about  the  mean.  Therefore,  unanswered  questions  include  how  close  to  an 
injection  slot  these  relations  are  useful,  and  where  and  how  the  concentration  fluc¬ 
tuations  must  be  considered. 


BOUNDARY-LAYER  MIXING  IN  NEWTONIAN  FLUIDS 


In  the  two-dimensional,  flat-plate,  turbulent  boundary  layer,  any  quantity  such 
as  the  instantaneous  local  concentration  of  polymer,  c,  will  be  a  function  of  x  (dis¬ 
tance  from  the  leading  edge),  y  (height  above  the  wall),  z  (spanwise  position)  and 
time  t.  Thus  c  =  c(x,  y,  z,  t),  but  the  time  average  is  independent  of  z-i.e., 
c  =  c(x,y)-~when  the  injection  is  uniform  across  the  span. 

Let  m  be  the  steady  rate  of  mass  discharge  per  unit  span  of  polymer  solution, 
and  let  c,-  be  the  polymer  concentration  (with  c;-  and  c  in  ppm  by  weight)  in  the  dis¬ 
charged  solution.  For  all  concentrations  of  interest,  the  fluid  density  can  be  taken 
as  that  of  the  water,  p  (in  grams/cm3 ).  Thus,  if  m  is  in  gram/cm  sec,  the  polymer 
mass  discharge  rate  per  unit  span  is  10-6  mc;-  gram/cm  sec  and  the  solution  volume 
discharge  per  unit  span  is  Qt  -m/p  cm2/sec. 

If  the  x  component  of  the  instantaneous  velocity  at  a  point  is  u  =  u(x,y,  z,  t), 
the  conservation  equation  for  the  polymer  mass  is 


Qici  -  J0  u(x,  y,  z,  t)c(z,  y,  z,  t )  dy 

because  of  the  uniformity  of  fluid  density.  The  overbar  indicates  the  time  average 
value  which  must  be  independent  of  z.  Interchanging  the  order  of  averaging  and  in¬ 
tegration  is  clearly  permissible,  so  that 

QiCf  =  So  uc(x,y)dy  (1) 


3 


Fiedler  and  Head  (Ref.  3)  studied  the  far-downstream  turbulent  mixing  of  smoke 
injected  into  a  wind-tunnel  boundary  layer  and  found  that  the  smoke  was  essentially 
confined  to  the  turbulent  fluid  and  vice  versa.  The  same  description  is  expected  to 
be  valid  for  polymer  molecules  in  water,  because  molecular  diffusion  of  polymer 
molecules  across  the  intermittency  interface  is  expected  to  be  negligible.  Thus  if 
I(x,  y,  z ,  t)  is  a  turbulence  detector  signal  (generated,  for  example,  by  a  vorticity  probe) 
that  is  1  when  the  flow  at  a  point  is  turbulent  and  0  when  it  is  not,  then  c(x,  y,z  ,t) 
can  be  replaced  by  c(x,  y,  z ,  t)I(x,  y,  z,  t).  Hence,  uc  in  Eq.  1  can  be  replaced  by 
luc  =  yuc(  where  y  ^  7  is  the  turbulence  intermittency  function  and  uc 1  is  the 
turbulent-conditional  average  value  of  uc  at  a  point.  Equation  1  now  becomes 

QjC;  =  JQ  yuc1  dy  (2) 


The  similarity  laws  for  Tic1  have  not  been  studied  either  experimentally  or  theo¬ 
retically.  However,  it  can  be  expected  that  uct  will  be  nearly  equal  to  ul  cl  if  the 
fluctuations  of  u  and  c  in  the  turbulent  fluid  are  small  relative  to  their  respective 
means,  ul  and  cf .  Thus,  considering  only  the  turbulent  fluid,  let  u  =  ut  +  u',  where 
W  -  0  and  c  =  c*  +  c\  where  c'x  =  0  so  that 


uc{  =  utct  +  u'c'1  -u*  cl 


u’c '*  (t/'2r  c' 2r)‘/2 

(u'2*  c'^)/z  utct 


Even  if  u'  and  c'  were  “perfectly  correlated,”  i.e.,  if  \urc't\  =  (ur*t  c^l)/z ,  and  if  the 
relative  rms  amplitudes  of  u'  and  c'  were  each  about  0.25,  the  error  in  the  approxi¬ 
mation  uc 1  «  u‘  cl  would  be  only  about  6%.  Thus  the  approximation  seems  accept¬ 
able  and  with  it,  Eq.  2  gives  (since  yc 1  =  c) 

Qici =  So  u‘ d  dy  (3) 


The  flows  considered  are  now  further  restricted  in  order  to  make  use  of  the 
similarity  laws  of  “equilibrium  boundary  layers”  (Ref.  4,  7)  for  large  Reynolds 
numbers.  This  makes  it  possible  to  approximate  the  u  distribution  throughout  the 
boundary  layer  by  using  the  velocity-defect  similarity  law  of  Coles  (Ref.  8): 


( U-u)/uT  =-(l/fc)lnT?  +  (n/fc)[2-vv(i?)]  (4) 

where  n  is  a  constant  in  a  given  type  of  equilibrium  boundary  layer;  17  ^  y/5n ;  Sn  is 
the  definition  of  boundary-layer  thickness  associated  with  Coles’  law  (with  subscript 
i r  to  distinguish  it  from  5  corresponding  to  y  for  u/U  =  0.99);  k  is  the  von  Karman 
“universal”  constant,  about  0.41,  and  believed  to  be  unchanged  with  friction  reduc¬ 
tion  caused  by  dilute  polymer  solutions;  and  w(7?)  is  given  by  a  tabulated  function 
(Ref.  8)  or  by  2  sin2  (7Tt?/2). 
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In  particular,  the  work  is  restricted  to  the  constant-pressure  type  of  equilibrium 
boundary  layer,  corresponding  to  that  on  a  flat  plate  at  zero  incidence  to  a  uniform 
free  stream.  It  is  thus  possible  to  make  use  of  recent  experimental  measurements 
which  suggest  the  similarity  laws  for  ul  and  c . 

Kibens  (Ref.  9)  measured  ul  in  a  nearly  constant-pressure  flow.  His  results  for 
u1 ,  u,  and  7  are  reproduced  in  Fig.  2,  where  it  can  be  seen  that  u{  is  at  most  about 
2%  less  than  u  and  then  only  where  7  (and  hence  c)  is  nearly  zero.  This  result  is 
quite  different  from  the  suggestion  of  Sarnecki  (Ref.  10)  that  u 1  should  satisfy  the 
logarithmic  law  of  the  wall  for  u  (cf.  Fig.  2).  With  the  approximation  u1  -u  in  Eq.  3, 

Qici  =  /cf  uc  dy  (5) 

Poreh  and  Cermak  (Ref.  1 1)  studied  the  two-dimensional  turbulent  mixing  of 
ammonia  gas  from  a  wall  line  source  under  a  wind-tunnel  boundary  layer.  They 
found  four  zones  of  development  of  the  mean  concentration  profile,  as  sketched  in 
Fig.  3.  In  the  fourth  and  final  zone,  they  found  that  the  normalized  distribution 
c(x,  y)/cw(x),  within  the  experimental  scatter,  was  a  function  of y/X  alone,  where 
A  is  defined  such  that  c/cw  =  0.5  at  y  =  A.  They  also  concluded  that  A/5  =  0.64  in 
the  final  zone.  Based  on  an  eddy-diffusion  analysis,  Morkovin  (Ref.  12)  described 
their  collected  data  for  c/cw  versus  y/X  in  that  final  zone  with  an  expression  conve¬ 
nient  for  the  purposes  of  this  report: 


c/cw  =exp[-0.693O/A)215]  (6) 

For  the  U/u7  range  of  the  experiments  of  Ref.  1 1,  it  can  be  shown1  that  5  «  0.95^, 
so  that,  assuming  that  A/5w  is  constant  and  independent  of  U/uT,  X/8n  0.58  is  the 
final  zone  law,  and  Eq.  6  becomes 


c/cw  =  exp[-2.270/6„)2-15]  (7) 

In  order  to  use  Eq.  4  and  7  in  Eq.  5,  the  latter  can  be  rewritten  as 

QiC i  =  uTcw  8n  /0°°  [(J/uT  -{U- u)/uT]c/cw  dr} 
or 

U8ncw/Qici  =  [/Q  c/cwdr\  -  (uT/U)  f*  (c/cw  )(U  -  u)/uT  cfy] 


^he  Gj  wake  function  of  Clauser  (Ref.  7)  leads  to  a  better  description  of  theu  profile  near  u/U  =  0.99 
than  does  vv(rj).  Using  Eq.  4  and  the  equivalent  expression  in  terms  of  Gj(v/6q),  where  6q  is  the  boundary- 
layer  thickness  used  by  Clauser,  it  can  be  shown  that  5q/6w=  1.08  if  fc  =  0.41.  Thus,  6/6^  =  1.08  6/5q  -  and 
6/8 Cl  can  1,0  dstenmned  as  a  function  of  U/uT  using  the  velocity-defect  law  plot  for  C \(y/f>Q\)- 
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Using  Eq.  4  and  7  to  evaluate  the  integrals,  while  taking  k  -  0.41  and  II  =  0.55  (the 
value  suggested  for  constant-pressure  flow  in  Ref.  8),  the  result  is 

U87!cw/Qici  =  (0.606  -  3.13ut/U)'1  (8) 

Before  comparing  this  relation  with  experiment,  it  is  of  interest  to  compare  cw 
according  to  Eq.  8  with  two  underestimates.  First  consider  ct ,  the  boundary -layer 
volume-flux  average  value  of  c,  which  disregards  the  concentration  profile,  so  that 

QjC'  =  Qtc,  (9) 

Since  only  the  turbulent  fluid  has  nonzero  c,  it  is  appropriate  to  make  Qt  the  turbu¬ 
lent  fluid  volume  flux  through  the  boundary  layer: 

Qt  k  ul  dy  =  /0“  yul  dy 

Using  the  approximation  ul  -  u  in  the  Qt  integral,  Eq.  9  gives 


Qici  =  UbnCt  [/(f  ydrl  -  iuTlU)  So  (10) 

It  is  shown  in  the  Appendix  that  the  assumption  that  y  =  y(rj)  with  the  experimental 
estimate  of  7(17)  for  constant-pressure  flow  from  Ref.  8  leads  to 

U8nct/Qici  =  (0.825  -  3.68 uT/U)A  (1 1) 

Thus  ct  «  0.1  cw  for  typical  values  of  U/uT,  so  that  c,  as  an  estimate  of  cw  would  be 
typically  about  30%  low. 

The  second  and  more  crude  estimate  for  cw  is  Cq  ,  which  is  obtained  by  replac¬ 
ing  Qt  in  Eq.  9  by  the  common  measure  of  volume  flux  through  the  boundary  layer,  Q 


Q=  So  u(dy/du)du 


This  estimate  neglects  both  the  concentration  profile  and  the  intermittency.  Using 
the  u  profile  law  assumed  in  this  work,  the  result  is 

U8ncQ/QiCi  =  (\  -3.78^/t/)"1 

which  yields  Cq  ~  0.6cw  for  the  high  U/uT  values  of  interest,  so  that  cQ  is  about 
40%  low. 

Figure  4  shows  the  plots  of  Eq.  8  and  1 1  versus  U/uT,  along  with  the  experi¬ 
mental  data  of  Ref.  1 1  for  U8cw/Qici.  Since  the  values  of  U/uT  for  the  experimental 
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FIG.  4.  Predicted  Mixing  Similarity  Laws  for  cw  and  ct  Versus  U/uT  (Eq.  8  and  11), 
Compared  With  Experimental  Data  for  U8cw/QjCj  (Ref.  11). 


data  are  not  given  in  Ref.  1 1,  they  were  estimated  by  means  of  the  relation  of  Schulz- 
Grunow  (Ref.  13)  for  the  experimental  U,  x,  and  v  values.  As  previously  mentioned, 

5  «  0.95w  in  the  experimental  U/uT  range,  so  that  the  center  of  gravity  of  the  experi¬ 
mental  points  would  be  roughly  on  the  theoretical  line  for  U8ncw/Qici  if  that  correc¬ 
tion  were  made.  A  best-fit  straight  line  through  the  data  would  have  a  steeper  slope 
than  the  theoretical  curve,  but  the  experimental  scatter  makes  this  conflict  of  little 
certainty  at  this  time.  It  is  interesting  that  Poreh  and  Cermak  fitted  the  line  of  zero 
slope  shown  in  Fig.  4  in  keeping  with  their  assumption  of  Eq.  5  with  a  power-law 
profile  u/U  =  (y/S)1/",  which  implies  that  U8cw/Qjci  is  constant. 

Equations  7  and  8,  for  c/cw  and  US^c^jQjCj,  seem  to  be  the  best  available 
similarity  estimates  for  the  far-downstream  mixing  of  a  passive  scalar  injected  uni¬ 
formly  from  a  line  source  into  a  two-dimensional,  constant-pressure  boundary  layer. 
The  generalization  of  these  results  to  include  the  case  of  polymeric  friction  reduction 
follows. 


THE  “NEGATIVE  ROUGHNESS”  ANALOGY 

Wall  roughness  does  not  enter  into  the  work  discussed  above,  since  wall  rough¬ 
ness  affects  the  skin  friction  but  not  the  velocity-defect  similarity  law  nor  the  inter- 
mittency  distribution,  y(y/8),  according  to  Ref.  14.  There  is  also  the  expectation2 


2  Private  communication  from  J.  E.  Cermak,  University  of  Colorado. 
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that  wall  roughness  does  not  affect  the  final  zone  law  for  c/cw .  It  is  notable,  there¬ 
fore,  that  polymeric  friction  reduction  has  been  likened  to  a  negative  roughness  on  the 
basis  of  the  following  comparison  of  characteristics: 

Wall  Roughness  Polymeric  Friction  Reduction 
tw  change:  >0  <0 

A B  change:  <0  >0 

where  tw  is  the  wall  shear  stress  pu 3,  and  A B  is  the  shift  of  the  constant  B  in  the  log¬ 
arithmic  law  of  the  wall: 

u/uT  =  (1  /k)  In (yuT/v)  +  B;  B  =  B0  +  A B  (12) 

where  B0  is  the  smooth-wall  value  for  Newtonian  fluids.  Since  it  is  known  that  the 
polymer  can  affect  the  flow  away  from  the  wall  (Ref.  15),  it  is  clear  that  the  analogy 
can  be  only  a  partial  one.  It  can  be  predicted  that  in  many  practical  cases  the  lower 
concentrations  away  from  the  wall  and  the  lower  Reynolds  stresses  there  render  the 
viscoelastic  effects  of  the  polymer  unimportant.  For  such  cases,  it  is  reasonable  to 
expect  that  the  velocity-defect  law,  the  intermittency  distribution,  and  the  mixing 
law  will  be  unchanged  with  polymeric  friction  reduction  since  they  are  unchanged 
with  wall  roughness.  The  similarity  relations  for  c/cw,  X/5^,  and  U8ncw/Qicj  should 
therefore  be  unchanged  with  polymeric  friction  reduction.  A  test  of  this  expectation 
is  made  below. 


A  COMPARISON  WITH  EXPERIMENT 


A  project  at  the  St.  Anthony  Falls  Hydraulic  Laboratory  (SAFHL)  of  the 
University  of  Minnesota  involves  the  mixing  process  in  a  turbulent  boundary  layer 
with  polymeric  friction  reduction.  The  investigation  is  being  made  in  a  9-foot-wide, 
open-water  channel  in  which  there  is  a  steady,  approximately  two-dimensional  flow 
of  about  1  foot  in  depth.  In  order  to  study  mixing  with  and  without  friction  reduc¬ 
tion,  dyed  polymer  solution  (Polyox  WSR-301)  or  dyed  water  is  injected  through  a 
tangential  wall  slot  extending  across  the  channel  near  the  effective  start  of  the  bound¬ 
ary  layer.  The  mean  dye  concentration  at  various  points  in  the  flow  downstream  of 
the  slot  is  determined  by  sampling  the  flow  through  pitot  tubes  that  are  also  used  to 
determine  the  mean  velocity  profile.  Two  representative  profiles  of  U  and  dye  c  for 
injection  of  water  and  500-ppm  Polyox  WSR-301  solution3  are  shown  in  Fig.  5  and 
6.  In  both  cases,  the  sampling  station  is  16  feet  downstream  of  the  injection  slot. 

The  values  of  Qt  and  dye  ct  are  listed  in  the  figures.  It  can  be  shown  that  16  feet  is 
far  enough  for  “final-zone”  behavior  to  be  expected  according  to  the  results  of 
Ref.  1 1 .  Since  the  dye  dilution  factor  at  the  wall,  ci/cw ,  is  about  5/0.06  in  the  case 


3  Private  communication  from  Prof.  J.  M.  Wetzel,  SAFHL. 
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of  polymer-solution  injection,  the  mean  polymer-solution  concentration  at  the  wall 
is  about  (0.06/5)  500  ppm  or  only  6  ppm.  This  low  value  and  the  moderate  wall 
shear  stress  (discussed  later)  suggest  that  the  negative  roughness  analogy  should  apply. 

Figure  7  shows  the  plots  of  y/X  versus  c/cw ,  based  on  Fig.  5  and  6,  plus  the 
prediction  of  Eq.  6.  The  comparison  is  quite  good  and  supports  the  expectation  of 
final  zone  behavior  and  the  negative-roughness  analogy. 


I  IG.  7.  Profiles  of  v/\  Versus  c/cw  from  Fig.  5 
and  6,  Compared  With  the  Final  Zone  Kmpirical 
Law. 


In  order  to  check  on  the  predictions  of  cw  and  X,  wT,and  5^  or  5  are  needed.  The 
values  of  uT  listed  in  Fig.  5  and  6  are  SAFHL  estimates  that  were  inferred  from  the 
velocity  profiles  by  use  of  the  logarithmic  law  of  the  wall.  These  rather  uncertain 
estimates  were  used  in  comparing  the  experimental  results  for  cw  and  X  with  theory. 

In  view  of  the  limited  amount  of  u  data,  5^  was  assigned  as  follows.  Each  displace¬ 
ment  thickness,  5*,  was  first  determined  from  the  corresponding  U  profile,  using 

U8*=  So  ( U-u)dy 


(13) 


The  SAFHL  data  fairings  shown  in  Fig.  5  and  6  were  used  in  Eq.  13  and  were  extrap¬ 
olated  to  u  =  0  using  the  logarithmic  law  of  the  wall.  The  resultant  5*  estimates,  the 
corresponding  values  of  U8*/v  (based  on  water  viscosity  at  the  test  temperatures),  and 
the  SAFHL  estimate  of  U/uT  for  each  injection  case  are: 


Water 

Polymer  Solution 

8*,  ft 

0.0389 

0.0287 

Uh*!v 

68,500 

26,900 

U/uT 

27.9 

30.5 

The  5^  values  were  calculated  from  the  5*  values  by  means  of  the  relation  for  8*/8n 
from  the  velocity-defect  law  of  Eq.  4  (for  k  -  0.41  and  II  -  0.55)  and  Eq.  13: 

8*/8n  =  (wT/C/)( n  +  1)1  k  =  3.78mt/C/  (14) 

Before  using  the  values,  a  check  on  the  consistency  of  the  t//wT  and  U8*lv 
values  can  be  made  by  use  of  the  relation  between  U/uT  and  U8*/v  for  constant- 
pressure  flow  of  Newtonian  fluids.  Figure  8  (reproduced  in  most  part  from  Ref.  4) 
gives  the  plot  of  U/uT  versus  \og(U8*/v)  for  Newtonian  constant-pressure  flow  on  a 
smooth  wall.  The  experimental  (aerodynamic)  data  of  two  investigations  and  two 


log  (U8*jv) 
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IMG.  8.  Comparison  of  Newtonian-Fluid  Relation  for  U/uT  Versus  log(t/S*/n) 
for  Constant-Pressure  Flow  Over  a  Smooth  Flat  Plate  (from  Ref.  4)  With  Data 
Points  Based  on  u  Plots  in  Fig.  5  and  6. 


semi-theoretical  expressions  are  shown,  along  with  the  two  SAFHL  data  points.  The 
solid  line  (from  Ref.  4)  corresponds  to  U/uT  =  5.75  log(U8*/v)  +  3.7.  The  dashed 
line  is  obtained  from  the  matching  of  the  velocity-defect  law  of  Eq.  4  with  the  log¬ 
arithmic  law  of  the  wall  (Eq.  12),  giving 

U/uT  =  ( 1  Ik)  In  (uT  8.  M  +  B  +  m/k 
By  use  of  Eq.  14,  this  becomes 

U/uT  =(l/k)\n(US*/v)  +  B  +  2Il/k-(\/k)\n[(n+  1  )/k]  (15) 

The  dashed  line  in  Fig.  8  is  U/uT  =  5.62  log((/5*/n)  +  4.44  according  to  Eq.  15  for 
k  =  0.41,  n  =  0.55,  and  the  Newtonian  fluid,  smooth-wall  value  of  B  ( B0  =  5.0). 

Both  of  the  semi-theoretical  lines  are  in  reasonable  agreement  with  the  aerodynamic 
data,  but  a  lower  slope,  i.e.,  a  smaller  value  of  (1  /k)  In  10,  may  be  needed  at  the  high 
end  of  the  U8*/v  range. 

The  fact  that  both  SAFHL  points  are  outside  the  spread  of  the  aerodynamic 
data  is  not  surprising.  The  low  U/uT  value  for  the  water-injection  case  (shown  by  the 
circle)  is  reasonable  because  a  A B  of  -3  corresponds  to  an  equivalent  sand-grain 
roughness  for  the  channel  bottom  of  a  likely  magnitude,  according  to  the  Prandtl- 
Schlichting  relation  (Ref.  16).  However,  an  overestimate  of  uT  could  produce  the 
same  result  for  a  smooth  wall.  The  high  U/uT  for  the  polymer-solution  injection  case 
(shown  by  the  triangle)  implies  a  A B  due  to  the  polymer  of  +5  if  the  reference  fric¬ 
tion  line  passes  through  the  water-injection  point,  parallel  to  the  smooth-wall  line. 
This  A B  is  very  low  in  comparison  with  the  pipe-flow  A B  values  found  by  Goren  and 
Norbury  (Ref.  17)  for  the  same  polymer  type  and  concentration.  If  the  wall  shear 
stress,  pw  2 ,  for  the  polymer-solution  case  is  estimated  as  300  dyne/cm2 ,  then  the 
extrapolation  of  the  A B  versus  tw  data  in  Ref.  17  according  to  the  Meyer  Law, 

A B  =  (a/2)  ln(rw/rM,  crit),  implies  a  A B  of  about  25  for  6  ppm  of  Polyox  WSR-301. 
The  lower  A B  may  reflect  any  of  several  factors,  such  as  wall  roughness,  an  over¬ 
estimate  of  uT,  a  breakdown  of  Meyer’s  relation,  or  differences  of  solution  properties. 

Bearing  in  mind  the  uncertainty  in  the  uT  and  8n  estimates,  the  experimental 
values  of  cw  and  5  can  be  compared  with  the  predictions  of  Eq.  8  and  X/ 8n  =  0.58. 
The  experimental  profiles  in  Fig.  5  and  6  and  the  previously  listed  5*  values  yield: 


Water 

Polymer  Solution 

U&*cw/Q'ci 

0.28 

0.26 

X,  ft 

0.23 

0.18 

The  values  of  8*/8n  according  to  Eq.  14  and  the  SAFHL  estimates  of  U/uT  were 
used  to  convert  these  to  values  of  U8ncw/Qjcj  and  X/5^.  The  results  are  as  follows: 
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Water 

Polymer  Solution 

0.136 

0.124 

UbjjQiCj 

2.07 

2.10 

0.81 

0.77 

These  values  are  shown  in  Fig.  9  along  with  the  predictions  of  the  scaling  relations. 

The  comparison  in  Fig.  9  supports  the  negative-roughness  analogy  in  that  the 
results  for  water  and  polymer-solution  injection  are  again  quite  similar.  However, 
while  the  cw  comparison  between  prediction  and  experiment  is  good,  the  large  per¬ 
centage  of  disagreement  for  X  shows  that  the  cw  agreement  is  partly  accidental,  since 
the  lower  predicted  X  values  mean  a  lower  predicted  flux  through  the  boundary  layer. 
This  conflict  reflects  the  fact  that  the  experimental  U  and  c  plots  in  Fig.  5  and  6 
yield  values  of  U  c  dy  which  are  50  to  60%  higher  than  the  experimental  values 
of  QjC^.  It  would  seem  that  either  the  flow  in  the  SAFHL  channel  is  significantly 
three-dimensional  or  that  the  approximation  of  uc  by  u  c  in  the  theory  is  in  signifi¬ 
cant  error.  The  latter  case  would  be  the  more  surprising,  since  the  difficulties  of  ob¬ 
taining  two-dimensional  flow  are  well  known,  and,  in  fact,  flow  visualization  in  the 
SAFHL  channel  showed  that  there  is  an  appreciable  lateral  contraction  of  the  dyed 
fluid  from  the  slot. 

In  general,  it  is  more  desirable  to  work  with  displacement  or  momentum  thick¬ 
ness  (5*  or  0)  instead  of  with  8n .  The  predicted  scaling  relations  in  terms  of  8*,  from 
Eq.  8  and  14,  are 


U8*cw/QjCj  =  3.78  (0.60 6U/uT  -  3.1 3)"1  (16) 

X/5*  =  0.\5U/ut  (17) 

Since  0  is  directly  involved  in  friction  reduction  work  and  since  its  experimental  de¬ 
termination  is  less  sensitive  to  error  in  u  near  the  wall,  it  is  the  most  desirable  thick¬ 
ness  parameter.  For  k  =  0.41  and  n  =  0.55,  use  of  the  velocity-defect  law  in  Eq.  4 
across  the  entire  boundary  layer  yields  [from  0  &  ( l/(/)2/J°  u(U -  u) dy] 

0/5*=  1-6.64  uT/U  (18) 

Thus  the  predicted  scaling  relations  in  terms  of  0  are 

UOcJQfCt  =  3.78  (0.606U/ut  -  3.13)"]  (1  -6.64  uT/U)  (19) 

X/0  =  0.15(//wr(l  -  6MutIU)~x  (20) 

As  an  example  of  the  use  of  these  predictions.  Fig.  10  gives  the  plots  of  c  versus 
y  according  to  Eq.  19  and  20  for  the  experimental  values  of  U ,  Q{,  ch  Ulur,  and  0 
from  Fig.  5  and  6.  The  0  values  determined  from  the  experimental  u  plots  are  listed 
in  Fig.  10;  the  dashed  curves  are  the  experimental  fairings  in  Fig.  5  and  6. 
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FIG.  9.  Comparison  of  Predicted  Law  for  US7!'cw/QjCi  and  \/6n  Versus  U/uT 
With  Data  Based  on  c  Plots  in  Fig.  5  and  6. 
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FIG.  10.  Comparison  of  Predicted  and  F.xperimental  c  Profiles  in  the  SAFHL  Tests, 
Based  on  Experimental  6  Values. 


This  dimensional  comparison  of  prediction  with  experiment  shows  that  there  is 
qualitative  agreement  of  the  c  profiles  but  a  sizable  underestimate  of  c  all  across  the 
profile.  This  is  due,  at  least  partly,  to  the  large  experimental  difference  in  the  SAFHL 
results  between  Qici  and  the  u  c  integral,  which  should  be  very  nearly  equal,  as  far  as 
is  known.  Thus,  if  most  of  the  difference  is  due  to  three-dimensionality,  the  quanti¬ 
tative  agreement  is  expected  to  be  better  for  truly  two-dimensional  flows. 

More  extensive  comparisons  of  the  predictions  of  Eq.  19  and  20  with  additional 
preliminary  data  from  SAFHL4  are  given  in  Fig.  1 1  and  12.  In  Fig.  1 1,  eleven  experi¬ 
mental  values  of  cw  for  x  =  16  or  40  feet  from  the  injection  slot  are  used  in  two  ways. 
The  unflagged  symbols  represent  experimental  values  of  (9,  ;  the  flagged  symbols  repre¬ 
sent  values  of  Qt  computed  from  Eq.  5,  using  the  experimental  u  and  c  distributions.5 
This  estimated  correction  for  the  effects  of  three-dimensionality  makes  the  experi¬ 
mental  points  fall  roughly  as  much  too  low  as  they  fell  too  high  before  the  correction; 
thus,  the  correction  is  not  adequate.  The  comparison  for  A/0  in  Fig.  1 2  is  not  affected 
by  the  Qt-  correction.  It  shows,  as  was  the  case  in  Fig.  9,  that  the  experimental  X  val¬ 
ues  arc  significantly  high,  i.e.,  by  about  30%  for  their  center  of  mass.  In  view  of  the 
small  number  and  the  scatter  of  the  data,  no  significant  differences  in  the  compari¬ 
sons  for  water  versus  polymer  solution  injection,  or  for  x  =  40  versus  1 6  feet,  can  be 
inferred  from  Fig.  1 1  and  12,  despite  the  two  exceptional  points  in  Fig.  12. 

The  amount  of  disagreement  between  the  theoretical  and  the  preliminary  ex¬ 
perimental  results  cannot  be  explained  at  this  time.  There  are  many  probable  causes, 
but  only  those  on  the  theoretical  side  are  discussed  here.  First,  the  empirical  simi¬ 
larity  laws  assumed  in  the  theory  are  subject  to  future  refinement.  For  example,  the 
A/5  =  0.64  value  of  Poreh  and  Cermak  for  the  final  zone  may  be  somewhat  low,  as 
consideration  of  Fig.  7  in  Ref.  1 1  seems  to  suggest.  Increasing  A/5  (and  hence  A/5^) 
would  tend  to  improve  the  agreement  in  both  Fig.  1 1  and  12,  since  a  higher  A  means 
that  a  lower  cw  is  required  by  the  polymer  conservation  equation. 

It  is  of  interest,  therefore,  to  see  what  shift  in  the  theoretical  relations  results 
from  a  conceivable  increase  in  A/5,  e.g.,  from  0.64  to  0.70.  It  can  be  seen  that  it 
causes  a  9%  increase  in  A/0  and  a  decrease  in  U6cw/QjCj  of  about  the  same  amount. 
These  shifts  would  move  the  theoretical  curves  noticeably  toward  the  centers  of 
mass  of  the  experimental  points  in  Fig.  1 1  (using  the  flagged  points,  of  course)  and 
in  Fig.  1 2.  Thus,  the  possibility  of  an  improved  estimate  for  X/8n  in  the  final  zone 
should  be  kept  in  mind. 

An  increase  in  X/8n  would,  however,  aggravate  a  questionable  feature  in  the 
final  zone  profile  of  c'/cw  which  is  implied  in,  but  not  of  direct  importance  to,  the 
present  analysis.  Because  the  scalar  concentration  in  the  nonturbulent  fluid  is  every¬ 
where  zero,  c,/cw  must  equal  ( c/cw  )/y.  Since  both  c/cw  and  y  approach  zero  at  the 


4 1  rom  Dr.  Wetzel;  sec  footnote  3,  page  9. 

5  Justin  McCarthy  of  Naval  Ship  Research  and  Development  Center  supplied  the  integrations  of  u  c. 
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edge  of  the  boundary  layer,  their  ratio  there  is  sensitive  to  error  in  either  part.  In 
particular,  the  c/cw  and  7  distributions  used  in  this  work,  which  were  taken  from 
separate  experimental  studies,  imply  that  cl/cw  dips  to  about  0.4  at  yl ~  0.75  and 
then  climbs  to  about  1 .0  at  y/8n  -  1 .0.  This  reversal  seems  peculiar,  and  the  con¬ 
ceivable  increase  in  X/S  from  0.64  to  0.70,  as  discussed  above,  would  by  itself  make 
the  reversal  more  pronounced.  Perhaps,  however,  the  cause  of  the  reversal  is  primarily 
an  overestimate  of  c/cw  for  larger  y/X  by  the  assumed  c/cw  profile.  In  fact,  inspec¬ 
tion  of  Fig.  1(a)  in  Ref.  12  suggests  that  such  is  the  case. 

Thus,  possible  improvements  in  both  final-zone  similarity  laws,  the  value  of 
X/5^  and  the  profile  of  c/cw  versus  y/X,  are  of  interest.  However,  since  small  changes 
in  c/cw  for  y/X  >  1.2  will  have  little  effect  on  the  present  prediction,  because  of  the 
small  contribution  to  the  u  c  integral  in  that  range,  the  possible  correction  to  X/6^  is 
more  likely  to  be  of  importance. 


SUMMARY 


The  available  experimental  data  and  empirical  relations  for  the  far-downstream 
mixing  of  a  passive  scalar  injected  into  a  Newtonian-fluid,  constant-pressure,  flat- 
plate,  turbulent  boundary  layer  appear  to  be  consistent  with  the  outer-layer  simi¬ 
larity  law  of  that  class  of  equilibrium  flows.  Thus,  invoking  the  analogy  of  polymeric 
friction  reduction  to  negative  roughness,  so  that  the  outer-layer  mean-velocity  simi¬ 
larity  (and,  presumably,  any  other  important  feature  of  the  outer  layer)  is  unaffected 
by  friction  reduction,  the  similarity  law  of  mixing  with  polymeric  friction  reduction 
is  predicted  to  be  the  same  as  it  is  without  polymeric  friction  reduction. 

In  order  to  develop  the  similarity  laws  within  the  framework  of  equilibrium 
boundary-layer  theory,  the  velocity-defect  law  of  Coles  is  adopted,  along  with  his 
definition  of  boundary  thickness  (5^ )  and  his  values  for  the  law-of-the-wak  •  param¬ 
eter  n  of  0.55  for  the  case  of  constant-pressure  flow. 

The  empirical  law  according  to  Poreh  and  Cermak,  as  formulated  by  Morkovin, 
of  the  normalized  mean  concentration  profile  recast  as  c/cw  versus  y/5n  for  the  final 
mixing  zone  is  adopted.  Their  scaling  law  for  wall  concentration  (cw )  is  replaced  by 
one  based  on  the  approximation  uc  ^Uc  (partly  justified  by  the  work  of  Kibens  and 
of  Fiedler  and  Head)  and  the  velocity  defect  law  of  Coles.  The  result,  in  reasonable 
agreement  with  the  data  of  Poreh  and  Cermak,  gives  US^/QjCj,  Ub*cw/Qici,  and 
U0cw/Qjci  as  slowly  varying  functions  of  U/uT.  A  noteworthy  prediction  is  that  cQ  , 
the  ordinary  boundary -layer  volume-flux  average  concentration,  is  about  0.6cw, ,  so 
that  Cq  ,  as  an  estimate  of  cw ,  is  still  about  40%  low. 

A  comparison  with  preliminary  data  on  mixing  in  an  open-channel  boundary 
layer  injected  with  dyed  water  or  dyed  Polyox  WSR-301  solution  suggests  that  the 
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predictions  of  the  negative  roughness  analogy  are  useful,  at  least  in  the  case  of  low 
polymer  concentration  and  moderate  wall  shear  stress.  Further  comparison  with 
more  nearly  two-dimensional  mixing  experiments  are  needed. 
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Appendix 


VOLUME-FLUX  AVERAGE  CONCENTRATION 


Following  upon  the  definition  of  ct  in  Eq.  9  in  terms  of  Qt,  Eq.  4  and  10  yield 


Qici/USnc,  =  fQ  ydr\  +  (uT/U)(\/k)  /J  y  In  r)drj 

+  (uT/U)(U/k)  /q  y(w  -  2)rfr?  (21) 

If  w  4  2  for  p  >  1 ,  then  Eq.  2 1  may  be  written  as 
QiCf/U^c,  =  [  1  -  (211  /k)ur/U]  Sq  7 dv  +  (ujUXl/k)  /„*  7  Inptfr? 

+  (uT/U)(U/k)  fQ  yw(r])dr}  (22) 


The  three  integrals  can  be  evaluated  if  7(17)  is  known.  There  is  a  sizable  amount  of 
experimental  information  on  7(.x,y)  in  Newtonian-fluid  boundary  layers.  Various 
investigators  have  found  that  their  distributions  can  be  described  accurately  in  terms 
of  two  parameters,  Y  and  o,  and  the  expression: 


1 


7  = 


exp 


dy 


=  (1/2)  erfc  I - —  j  [using  erf  (- a )  ^  -erf  (a)]  (23) 

\  / 

This  relation  has  a  simple,  but  not  unique,  interpretation:  if  for  any  given  station  x 
there  is  at  every  instant  a  single  interface  at  y  =  Y(x,  t)  between  turbulent  (y  <Y) 
and  irrotational  flow,  and  if  Y(x,  t )  has  a  Gaussian  probability  density  distribution 
with  mean  value  Y  and  standard  deviation  a,  then  the  turbulence  intermittency  dis¬ 
tribution  y(x,  y )  is  given  by  Eq.  23.  7£quals  0.5  at  y  =  Y  and  is  essentially  1  for 
y  <  Y  -  2.5 o  and  essentially  0  fory  >  Y  +  2.5 o.  Fiedler  and  Head  (Ref.  3)  found  in 
their  experiments  with  both  favorable  and  unfavorable  pressure  gradients  that6  Y /5 


6 The  meaning  of  6  in  Ref.  3  is  not  defined  but  was  assumed  to  be  the  99%  point  in  this  work. 
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and  0/6  appeared  to  depend  only  on  the  profile  form  factor,  H  &  8*16 ,  and  that 
there  was  no  systematic  effect  of  Reynolds  number  for  the  range  covered,  viz. 

U6/v*>  1000  to  4000. 

If  this  relation,  that  Y/8  and  a/5  depend  only  on  H,  were  exact,  so  that 
7  =  7(77,  H),  it  would  mean  that  the  intermittency  in  a  constant-pressure  boundary 
layer  along  a  line  of  constant  y/8  (ory/5ff,  to  be  precise)  would  vary  with  x  as  a 
resultof  the  variation  of  H  with  U/uT.  It  seems,  on  the  other  hand,  more  reasonable 
that  Y/5ot  and  a/5 v  should  be  as  nearly  independent  of  x  in  a  constant  pressure, 
equilibrium  flow  as  is  the  velocity-defect  (U -  u)/uT  versus  y/8n,  since  intermittency 
is  primarily  a  property  of  the  outer  layer. 

Several  investigators  have  measured  intermittency  distributions  in  nearly  constant- 
pressure  flow  and  obtained  very  similar  results.  In  particular,  Coles  (Ref.  8)  reanalyzed 
Klebanoff’s  data  (Ref.  18)  and  obtained  Y/8n  =  0.825  and  a/5 n  =  0.148,  after  assum¬ 
ing  II  =  0.55.  The  outer-layer  similarity  prediction  is  thus  thatJVS^  =  0.825  inde¬ 
pendent  of  U/ur.  Figure  13  contrasts  this  prediction  with  the  Y/8  prediction  from 


FIG.  13.  Two  Predictions  for  Y/6n  Versus  U/uT  for  a  Constant-Pressure 
Boundary  Layer. 


Ref.  3,  combined  with  the  well-known  relation  between  H  and  U/uT  for  large 
boundary-layer  Reynolds  numbers: 


21 


H  =  (l  -GuJU)'1 

G  4  Jq  [(U-  u)/ut  ] 2  dr) I /q  (U  -  u)!ut  dr\ 

G  =  6.64  (as  in  Eq.  18) 

(This  value  of  G  is  quite  reasonable,  as  can  be  seen  in  Ref.  4  and  14.) 

Since  the  difference  between  and  5  are  only  of  the  order  of  10%  for  ordinary 
U/uT  values,  the  disagreement  in  Fig.  13  cannot  be  attributed  to  that  difference.  The 
prediction  based  on  the  Y/8(H )  relation  of  Ref.  3  clearly  seems  more  suspect  than  the 
other  because  the  experiments  of  Fiedler  and  Head  for  low  H  values  were  for  very 
favorable  pressure  gradients.  Thus  the  constant-7/5^  behavior  is  assumed  here. 

Using  Y/8n  =  0.825  and  0/8^  =  0.148  and  Eq.  23,  the  integrals  in  Eq.  22  can  be 
evaluated.  The  first  integral  can  be  evaluated  exactly: 

f0~ldv  =  7(0,  Y/Sn,  a/8w)Y/8n  +  (o/8n^2t)  exp [-(f/o^/2)2] 

For  the  assumed  values  of  Y/8n  and  o/8v,  a  very  good  approximation  is 

f0  ydv  «  Y/Sn 

The  second  integral  was  evaluated  numerically  except  for  an  interval  containing  the 
singularity  in  lnrj,  which  was  handled  analytically.  The  result  is 

So  7  ln77cfy  «  -0.7727(0,  Y/8n,  o/8n)  -  0.198  «  -0.970 

Numerical  integration  of  the  third  integral  gives 


f0  yw(rj)dr]  0.668 

These  three  results  and  k  =  0.41  and  II  =  0.55  in  Eq.  22  yield  the  scaling  law 
for  ct  given  in  Eq.  1 1 . 
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